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Increasing attention has been recently devoted to protein-based nanobiosensors. The main reason is the huge number 
of possible technological applications, going from drug detection to cancer early diagnosis. Their operating model 
is based on the protein activation and the corresponding conformational change, due to the capture of an external 
molecule, the so-called ligand. Recent measurements, performed with different techniques on human 17-40 olfactory 
receptor, evidenced a very narrow window of response in respect of the odour concentration. This is a crucial point for 
understanding whether the use of this olfactory receptor as sensitive part of a nanobiosensor is a good choice. In this 
paper we investigate the topological and electrical properties of the human olfactory receptor 17-40 with the objective 
of providing a microscopic interpretation of available experiments. To this purpose, we model the protein by means of a 
graph able to capture the mean features of the 3D backbone structure. The graph is then associated with an equivalent 
impedance network, able to evaluate the impedance spectra of the olfactory receptor, in its native and activated state. 
We assume a topological origin of the different protein electrical responses to different ligand concentrations: In this 
perspective all the experimental data are collected and interpreted satisfactorily within a unified scheme, also useful 
for application to other proteins. 



?3 1 Introduction 



I 

o 
o 



> 

00 
O 



> 

X 



Olfaction is one of the most refined senses in mammals: 
It allows to recognize many thousands of different odours. 
In general, each odour is a mixture of different molecules, 
each of them able to activate a specific set of olfactory 
receptors (ORs), present in the nasal epithehum pfliniElll]. 
By contrast, each OR selectively recognizes few specific 
molecules (ligands) [5]. In such a way, the sequence of 
activated ORs identifies the odour univocally [7] . This 
mechanism indicates the ORs as the natural candidates for 
the production of a new generation of hybrid nanodevices 
for the detection of specific ligands 13 [IHl HH [HI [l3] . To 
this purpose, these nanobiosensors should integrate animal 
ORs with an electrical/electronic transducer [5|. 

Conventional electronic noses are biomimetic sensors 
able to detect specific odours. The detection is usually 
performed by means of solid-state gas sensors: They trans- 
late the captured information into a digital signal which 
is elaborated by appropriate algorithms. To date, the sen- 
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sitivity and selectivity of conventional devices are of good 
level but still very far from those of a mammalian nose, 
the gold standard of an olfactory biosensor [H [U] . The use 
of ORs as active part of a sensor aims to improve its selec- 
tivity, to reduce its size for achieving an on-field detection, 
to cut the costs and, as a future challenge, to obtain for 
each odour of interest a specific device [J [H [HI IH [HI HH US] • 

These preliminary remarks are mainly supported by the 
observation of a change of the electrical response in an 
OR-based devices [m [15l [16] . This change coincides with 
the protein activation due to the presence of the specific 
ligand. 

All the ORs are G-protein coupled receptors (GPCRs): 
This family of proteins includes receptors of many different 
molecules and also the light receptors (opsins). GPCRs 
have structure and function quite similar. Previous stud- 
ies on opsins [S] and rat OR 17 [H [T7] validated a the- 
oretical model able to describe the electrical properties 
of these proteins through an impedance network protein 
analogue (INPA). Recent investigations on human OR 17- 
40 [131 [T51 [12], valuable for a practical realization of a 
nanobiosensor, led to interesting results not yet fully un- 
derstood on a microscopic ground. In particular, measure- 
ments were performed with electrochemical impedance 
spectroscopy (EIS) on a two terminal device with the ac- 
tive part controlled by this OR: It was observed that the 



polarization resistance showed a peculiar bell-shaped be- 
haviour, at increasing values of the specific odorant con- 
centration [13]. The peak was centered at low concen- 
trations of the specific ligand (around 10""'^° M). Further- 
more, by using complementary techniques like the differen- 
tial surface plasmon resonance (SPR) and the differential 
bioluminescence response, in addition to the first peak, a 
second maximum response was observed at a higher con- 
centration of the odorant (around 10"^ M)[TH]. These re- 
sults were also confirmed by differential conductance mea- 
surements [16]. 

The aim of this paper is to group recent investigations 
performed on human OR 17-40, in particular those ob- 
tained from EIS measurements, into a unified microscopic 
model. By using the INPA approach we correlate the elec- 
trical response of the given protein to its topological prop- 
erties. Accordingly, it is possible to predict how the modi- 
fications of the latter, as induced by the capture of ligands, 
influence the former. The understanding of this correla- 
tion is promising for a future use of human OR 17-40 in 
nanobiosensors. Finally, we conjecture the existence of a 
correspondence between the collective response of an en- 
semble of proteins to a specific odour concentration, and 
the variation of the calculated single protein impedance 
at different values of the network free parameter Re. This 
parameter defines the maximal distance between two inter- 
acting amino acids. In this way, the single protein model 
is able to provide estimations on the electrical response of 
a protein based macroscopic device. 

The paper is organized as follows. Section [2] briefly sum- 
marizes the theoretical model. Sections [3] and |4] detail the 
main results of the paper and illustrate a conjecture for 
describing the protein modiflcations due to different lig- 
and concentrations. A conclusive section summarizes the 
major achievements and the perspectives of this research. 



2 Theory 

The INPA approach, as detailed in previous papers [3 [H] , 
correlates the electrical and topological properties of a 
sensing protein. The main objective of INPA is to pre- 
dict the protein electrical modifications as induced by a 
change of its topology. For the reader convenience, its 
main features are here briefly recalled. The approach is 
structured in two steps. 

The first step produces the protein analogue graph. The 
input data of the model are the tertiary (3D) structures of 
the protein. These structures are known from the protein 
data bank (PDB) available in literature [20] or similarly 
(see, for example, Ref.[17]). Given the 3D structure, the 
protein is represented by means of a topological network 



(a graph) . Each node of the graph corresponds to a single 
amino acid. If two amino acids are closer than an assigned 
interacting radius. Re, then a link is drawn between the 
corresponding nodes. 

The second step dresses the graph. By associating an 
elemental impedance with each link, the graph becomes 
an impedance network. In other words, each link mim- 
ics a privileged channel for charge transfer and/or charge 
polarization [T7] . 

The elementary impedance between nodes i and j, Zij, 
consists of an Ohmic resistance, i?, parallely connected 
with a parallel plate capacitor, C, filled with a dielectric 

[siiniiiz]: 
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where Ai.j = t^{Rc^ — IfJ^), is the cross-sectional area 
between two spheres of radius Re centered on the i-th and 
j-th node, respectively; Uj is the distance between these 
centers, p is the resistivity, taken to be the same for every 
amino acid, with the indicative value of p = lO'^" Jim; i = 
\/—l is the imaginary unit, eg is the vacuum permittivity, 
and io is the circular frequency of the applied voltage. The 
relative dielectric constant of the couple of i, j amino acids, 
e^ j, is expressed in terms of the intrinsic polarizability of 
the i,j amino acids. 

By positioning the input and output electrical contacts, 
respectively on the first and last node (corresponding to 
the first and last amino acid of the protein sequential 
structure) , the network is solved within a linear Kirchhoff 
scheme [9] [TOl [17]. Within this framework, the protein 
is mapped into an electrical network whose topology re- 
flects the protein conformation. When the conformation 
changes, the network topology and, in turn, the electrical 
response, follows this change. 



3 Material and Method 

To date, no crystallographic data are available for the 3D 
structures of the human OR 17-40. Accordingly, the 3D 
structures used in this paper are obtained from an homol- 
ogy modeling [IT] . Those are calculated by using the light 
receptor bovine rhodopsin in the native state (in dark), 
as template for the henceforth mentioned native human 
OR 17-40, and the bovine rhodopsin in the active state, 
metarhodopsin II (in light), as template for the henceforth 
mentioned active human OR 17-40. 

The INPA approach is very flexible, and allows one to 
carry out different kinds of investigations. In particular, 
here we point: 



i) To compare the contact maps of the protein, in the 
native and active state, as a function of the interacting 
radius. Re. These 2D maps visuahze the connected parts 
of the network, and are drawn by assigning to each couple 
of hnked amino acids, say i — j, a point of coordinates 
i,j. Two amino acids are hnked if their distance is less 
than i?cj otherwise they remain separated. Therefore, by 
comparing the contact maps of the protein in the native 
and active state we gain access to a direct image of the 
main changes in the protein topology. 

ii) To calculate the relative variation of resistance be- 
tween the native and active state of the OR. Since the 
elementary impedance depends on the distance between 
amino acids, when the protein conformation changes also 
the global impedance, and thus the resistance, changes. 

iii) To calculate, within the standard frequency range 
1 niHz -;- 100 kHz, the protein impedance spectra for the 
native and active states. In this way it is possible to draw 
the theoretical Nyquist plot that should be compared with 
the experimental one. 

Finally, the different electrical responses, as calculated 
from the networks pertaining to the native/active state, 
are here taken as an estimate of the measured electrical 
responses implied by the protein conformational change. 



4 Results 

The main results concerning sections [3] and [2] are detailed 
in the following subsections. 



4.1 Protein topology 

The global insight on the protein conformational change, 
as induced by the ligand capture, is given by the contact 
maps reported in Fig. [l] and Fig. [2J In these figures the 
diagonal is a symmetry axis that allows one to compare 
native and active states in the same figure, the native state 
is represented by points satisfying i < j while the active 
state by points i > j. Furthermore, the contact maps are 
calculated for two different interaction radii, specifically 
6 A in Fig. [ij and 20 A in Fig. [2| From these figures we 
can observe that, by increasing the value of Rc^ the density 
of points and thus the connectivity increase substantially, 
as expected. In particular, the main differences between 
the native and active state are in the increased connec- 
tivities among the helices h3-h4-h6, where the principal 
binding pocket is presumably located [2l] . Notably, signif- 
icant differences between native and active states survive 
at large R^ say up to Re values of about 50 A. 
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Figure 1: (Color on line) Contact maps of human OR 17- 
40 for Rc= 6 A. The data of the native state are reported 
on the left side of the diagonal (black full squares), the 
data of the active state are reported on the right side of 
the diagonal (green open squares). The ellipses signal the 
main differences between the native and active state which 
are attributed to the transmembrane helices h3-h4-h6. 

4.2 Protein resistance 

The value of the protein resistance in its native state is in 
general different from that in its active state, and the mag- 
nitude of this difference depends on Re as reported in Fig. 
3. Here the main result is the non-monotonic behaviour of 
the resistance change at increasing values of Re- The re- 
gion of maximum sensitivity to the conformational change 
is found for Re in the range 6 -^ 14 A. Furthermore, in the 
ranges of Re values 18 -;- 26 A and 38 -;- 65 A, calculations 
evidence an inversion of the resistance variation, with the 
active state becoming less resistive than the native state. 
In terms of the electrical network, such an inversion is in- 
terpreted as a stronger increase of parallel with respect to 
series connections of the active link resistances. 

Recent advances on the activation mechanism of these 
proteins [55] depict their dynamics in terms of a complex 
series of intermediate conformational transitions (that in- 
volve the disruption of non covalent intramolecular inter- 
actions) rather than in terms of a simple on/off switch. 
These interactions stabilize the protein state in a specific 
equilibrium condition, which depends on the ligand con- 
centration. An increasing concentration of ligands induces 
the disruption of stabilizing interactions, enabling the re- 
ceptor to evolve toward a more active state. 

In our model, we conjecture that the disruption of sta- 
bilizing interactions is related to a reduction of the in- 
teraction radius Re, which corresponds to a reduction of 
connections between amino acids. Accordingly, we assume 
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Figure 3: Relative resistance variation of human OR 17-40 
as function of the interacting radius Re ■ 



Figure 2: (Color on line) Contact maps of OR 17-40 for 
Rc= 20 A. The data of the native state are reported on 
the left side of the diagonal (black full squares) , the data 
of the active state are reported on the right side of the di- 
agonal (green open squares). The ellipses signal the main 
differences between the native and active state which are 
attributed to the transmembrane helices h3-h4-h6. 



the existence of a correlation between the interacting ra- 
dius, Re, and the concentration of a specific odorant, A. 
In this way, the change of the single protein resistance vs 
Re (see Fig. 3) is correlated with the analogue quantity 
measured for different ligand concentrations pJJ. In es- 
tablishing this correspondence we assume that the peak 
of the electrical response of human OR 17-40 (occurring 
at A = 10"^" M) corresponds to Rc~ 52 A, as resulting by 
Fig. 3. Furthermore, in this perspective. Fig. 3 gives the 
evolution of the relative resistance variation at different 
ligand concentrations. Starting from very large R^ where 
there is no difference between state responses (small A val- 
ues in our conjecture), at smaller Re the resolution of the 
electrical response evolves toward larger values (large A 
values). In doing so, the resolution exhibits at least two 
bumps. 

The functional dependence of Re on A is found to fol- 
low a power law, as evidenced by the dashed line in Fig. 
[3j To enforce the above conjecture. Fig. |4] reports Re as 
a function of a specific odorant concentration. Here, the 
data are obtained by comparing the experimental response 
for different helional concentrations [XT' with the theoret- 
ical relative resistance variation of the single protein at 
different Re values, as reported in Fig. |3] 

Figure [5] reports the polarization resistance variation as 
a function of the odorant concentration, as measured for 
the specific cases of heptanal and helional [14]. Theoret- 
ical results (dashed curve) are found to reproduce quali- 
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Figure 4: Interacting radius (Re) as function of the specific 
ligand helional concentration, A. Dashed line is a guide to 
the eye. 



tatively well the experiments (symbols) , which evidence a 
typical bell-shaped behaviour with a maximum sensitivity 
for the case of helional at a molar concentration of 10~^° 
M. By construction, the quantitative comparison between 
theory and experiments cannot discriminate the different 
sensitivity to different odorants of the OR. In any case, 
at this stage of investigation we believe that the overall 
agreement between theory and experiments is satisfactory 
enough to further support our conjecture. 

Finally, we notice that the dose-response depicted in 
Fig. 5 departs from the saturation behaviour observed 
for other receptors (like rat OR 17) where an increasing 
concentration of the specific ligand produces a systematic 
increasing of the protein response [T3] . According to our 
conjecture, the origin of this difference should be found 
in the peculiar modification of the topology undergone by 
human OR17-40. 



DC 

•^ 0.01 K 



-12 





1 


1 1 


' 


: 




♦ 






- 


r 


< 


t ' 


♦ 




^ 








i 


• 




!- 








^ 




theory 






'- 




• helional 






. 




n heptanal 


1 1 




" 



-10 



log A 



Figure 5: Per cent change of the polarization resistance 
at increasing concentration of specific odorants for human 
OR 17-40 at room temperature. Symbols refer to experi- 
ments, with bars giving the estimated uncertainty |151ll4j. 
dashed line refer to theoretical expectations. 



4.3 Protein impedance spectrum 

The impedance spectrum of the single protein is explored 
over a wide range of frequencies and the results are given 
by means of the associated Nyquist plot. This plot is ob- 
tained by drawing the negative imaginary part versus the 
real part of the global impedance, within a given frequency 
range (typically from 1 mHz to 100 kHz as in experiments 
[131 )• Figure [6] reports the Nyquist plots of human OR 
17-40 with the impedance normalized to the static value 
of the native state. Symbols pertain to experiments, with 
crosses referring to the absence of a specific odorant, full 
(empty) squares to the presence of the specific odorant he- 
lional (heptanal) at the concentration reported in the fig- 
ure. Curves pertain to theoretical results where the single 
protein is taken to be representative of the entire sample, 
and with continuous (dashed) lines referring to the native 
(active) state. Within our model, the Nyquist plots arc ob- 
tained by using as input data the networks corresponding 
to the native and active states at the cut-off radius which, 
according to Fig. [5] gives the maximum resolution. The 
agreement between theory and experiments is found to be 
satisfactory from a qualitative point of view and, apart for 
a significant underestimation of the maximum experimen- 
tal value, acceptable from a quantitative point of view. We 
remark that the near ideal semicircle shape of the exper- 
imental Nyquist plot is well reproduced, thus confirming 
that the network impedance model behaves closely to a 
single RC circuit as expected by the presence of a rather 
uniform distribution of time constants associated with the 
different values of the resistance and capacitance of the 



cc 



E 
N 




0.2 0.4 0.6 0.8 
ZRe/ZRe(0) 



Figure 6: Nyquist plot of human OR 17-40 in the ab- 
sence and in the presence of a specific ligand (heptanal, 
helional). The impedances are normalized to the static 
value of the native state, Z'^{Q) = 33 Kfi cm^. Symbols 
pertain to experiments with: Crosses referring to no odor- 
ant, empty (full) squares to an heptanal (helional) concen- 
tration of 10"^" M at room temperature [inilllj. Curves 
pertain to theoretical results with: Continuous curve re- 
ferring to the native state configuration as input data with 
Re = 70A, dashed line referring to the activate state con- 
figuration with Re = 46 A. 



links [17] . The scarcity of experimental data and the lack 
of a well certified knowledge for the 3D structures of the 
considered proteins lead us to consider these results as a 
first but significant step towards a microscopic modeling 
of the electrical properties of this olfactory receptor. 

5 Conclusions and perspectives 

Recent EIS experiments on human OR 17-40 showed that 
the mechanism of odorant capture can be monitored by 
means of the modification of the impedance spectra I 14j . 
These results are generally confirmed by surface plasmon 
resonance and bioluminescence response of appropriately 
prepared OR 17-40 samples [HldH]- The EIS results are 
here microscopically interpreted on the basis of the con- 
formational change of the protein tertiary structure, in- 
duced by the sensing action. Accordingly, the conforma- 
tional change induces a variation of the protein electrical 
response that is detectable with the EIS technique and 
can be reproduced by means of the INPA modeling. In 
particular, a possible correlation between the interaction 
radius, Re, at the basis of charge transfer between amino 
acids, and the odorant concentration, A, is inferred by 
comparing theory with experiments. 

From one side, the present INPA approach is rather es- 



sential and the qualitative agreement with available exper- 
iments is used as a proof of concept that is promising for 
the realization of an electronic nose based on an array of 
nanobiosensors. Furthermore, the narrow window of sen- 
sitivity of human OR 17-40 to the odorant concentration 
poses potential limitations in the advantages of receptor 
based biosensors in respect of traditional sensors. These 
limitations should be overcome by investigating similar 
ORs which are, like rat OR 17 [T7], potentially free from 
these drawbacks. 

From another side, relevant information on the protein 
tertiary structure and the mechanisms of charge transfer 
can be extracted from the microscopic interpretation of ex- 
periments. The next step is the production of new identi- 
fied ORs in efficient heterologous systems (mammalians or 
yeast cells), to evaluate their functional response. In par- 
allel, improvements of the theoretical model should take 
into account the role of temperature, hydration, and other 
details. 
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